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A simple model for nuclear structure functions in the region of small x and small and 
moderate Q^, is presented. It is a parameter-free extension, in the Glauber-Gribov approach 
to nuclear collisions, of a saturation model for the nucleon. A reasonable agreement with 
experimental data on ratios of nuclear structure functions is obtained. Nuclear effects in 
the longitudinal-to-transverse cross section ratios are found to be small. Predictions of the 
model for values of x smaller than those available to present experiments are given. The 
unintegrated gluon distribution and the behaviour of the saturation scale which result from 
this model are shown and discussed. 
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1 Introduction 



The differences between the structure functions measured in nucleons and nuclei 
the so-called EMC effect, are a very important feature for the study of nuclear structure 
and nuclear collisions. At small values of the Bjorken variable x (< 0.01, shadowing 
region), the structure function F2 per nucleon turns out to be smaller in nuclei than 
in a free nucleon. The nature of this shadowing is well understood qualitatively: In 
the rest frame of the nucleus, the incoming photon splits, at high enough Q^, into a 
qq pair long before reaching the nucleus, and this qq pair interacts with it with typical 
hadronic cross sections, which results in absorption ^, ^, |^, 0; in this way nuclear 
shadowing is a consequence of multiple scattering and is thus related with diffraction 
0. An equivalent explanation in the frame in which the nucleus is moving fast, is that 
gluon recombination due to the overlap of the gluon clouds from different nucleons, 
makes gluon density in nucleus with mass number A smaller than A times that in 
a free nucleon ||, ^. These studies have received great theoretical impulse with the 
development of semiclassical ideas in QCD and the appearance of non-linear equations 



for evolution in x in this framework ||T0|, |TT], 0, [T^, |T4|, |15| [T^, although saturation 
appears to be different from shadowing Jl^ (see [^] for a simple geometrical approach 
in this framework). 



On the other hand, a different approach is taken in |T^: parton densities inside 
the nucleus are parametrized at some scale ~ 2 GeV^ and then evolved using the 
DGLAP evolution equations. In this way, the origin of the differences of partons 
densities in nucleons with respect to nuclei is not addressed, but contained in the 
parametrization at Qq which is obtained from a fit to experimental data. 

The results from different models usually depend on additional semiphenomeno- 
logical assumptions and often contradict each other. For example, concerning the 
Q^-dependence of the effect, in ^ |^ it is argued that qq configurations of a large 
dimension give the dominant contribution to the absorption, which results essentially 
independent of (this is the case in until extremely small x, where a dependence 
cx In Q/Q appear related with the use of BFKL evolution [0). On the other hand, in 
the gluon recombination approach of the absorption is obtained as a clear higher- 
twist effect dying out at large Q^. Finally, in the models which use DGLAP, all 
Q^-dependence comes from QCD evolution and is thus of a logarithmic, leading-twist 
nature, see for a comparison between multiple scattering and DGLAP approaches. 
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Predictions (particularly for the gluon density) on the x-evolution towards small x turn 



out to be very different p3|, p4 |. 

In practice these studies are of uttermost importance to compute particle production 
in collisions involving nuclei. For example, in the framework of coUinear factorization 
25|| parton densities in the nucleus following the spirit of [|19| are needed, see e.g. ||2^ 



for recent applications to inclusive particle production in heavy ion collisions. While 
this scheme is suitable to compute particle production at scales Aqfj^ ^ ^ and 

thus in the hard region, for semihard (minijet) production, Aqj^^d -C s, the kx- 

factorization scheme should become the suitable one (in ^ applications 



to heavy ion collisions can be found). Here the tool that is needed is the so-called 



unintegrated gluon distribution (see a precise definition below and for a review of 
its situation for nucleons). Let us stress that different approaches give very different 
predictions for multiplicities at RHIC and LHC . 

In this paper we present a simple model for nuclear structure functions in the region 
of small X (< 0.01) and of small and moderate (< 20 GeV^) in the dipole picture 
P, It consists on an extension to nuclei, using the Glauber-Gribov picture ^4 



of the saturation model for the proton in without any new parameter (in 
a similar strategy is used but with a simplified dipole-nucleon cross section and no 
comparison with experimental data on nuclear structure functions is performed; also 
in [|l6l this ansatz is used but just as an initial condition for a non-linear evolution 



equation). The plan of the paper is as follows: In the next Section the model will be 
described. In Section ^ a comparison of the results of the model with available data 
on F2A will be shown. In Section ^ the unintegrated gluon density obtained in this 
model, together with the saturation scale it implies, will be discussed. Finally, in the 
last Section conclusions and possible applications of the model will be outlined. 



2 Description of the model 

The nuclear structure function F2 can be standardly defined via the cross sections (Jt,l 
for the collision of the transversal (T) or longitudinal (L) virtual photon of momentum 
q, = —Q^, on the nucleus A of momentum Ap: 

F2A{x,Q') = ^^{a^ + at); (1) 
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this expression holds for the small x region where we are going to work. Both cross 



sections can be conveniently presented in the dipole model 0, Q via the cross section 



o'dAiXjT), for the scattering of a colour dipole of transverse size r on the nucleus: 

(^tA^) = J d'^r pT,L{r)adA{,x,r), (2) 



where pt,l are the distributions of colour dipoles created by splitting of the incident 
photon into qq pairs 0] : 

□ 2 i 



PT(r) = '^T.4f + (1 - ")']«'K;(er) + mjKS(er)} (3) 

and 

PL(r) = ^g^E4 1'daa^il-afKlier). (4) 

Here summation goes over flavours (which will be limited to 3 or 4, see the comments 
below Eq. (|^)), = Q'^a{l — a) + m'j, and m/ and Zf are respectively the mass and 
electric charge in units of the proton charge e, of the quark of flavour /. 

For the total cross section of a dipole on a proton we will use the saturation model 
in p5[, which provides a good description of inclusive and diffractive experimental data 



on for x < 0.01 and small and moderate < 20 GeV^. Other models which also 
describe data in this region could be used, as that in |3^, see ||3^ for a comparison 



between this model and that of |3^, which we will employ because of its simplicity. 
The form of the cross section isQ 



adp[x,r) = (To 
with 



exp 



(5) 



Qli-) = Ql{^) , x = x(l + ^j. (6) 



In |^5| some parameters were kept fixed: Qq = 1 GeV^, m„ = = = 0.14 GeV 



and rric = 1.5 GeV. Three parameters were determined from a fit to data and have 
values (Jo = 23.03 (29.12) mb, A = 0.288 (0.277) and Xq = 3.04 ■ 10"^ (0.41 ■ 10"^) for 
the 3-flavour (4-flavour) version of the model. is called the saturation scale, which 
we will discuss in Section H. 



^In this model no impact parameter of the proton is cxphcitly given, so it cannot be considered 
a pure eikonahzation of some elementary amplitude |T7[ ; on the contrary, in our extension to nuclei, 
Eqs. M) and m) below, the impact parameter of the nucleus is explicitly taken into account. 
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Let us comment that, neglecting the influence of quark masses in variable e which 



has ||3^ small effect in the case of light quarks, this model implies, in its 3-flavour 
version, an exact scaling of the cross sections with r = Q"^ /Ql- This scaling has 
been shown to be fulfilled to a good approximation by all DIS data for x < 0.01 
(which go up to ~ 450 GeV^)0. So this model, which leads for large to Bjorken 
scaling, is apparently able to mimic along a wide Q^-region the QCD evolution (see 
39(1 for improvements of this model to include DGLAP evolution). This r-scaling has 



been argued to hold for such a large > (~ 1 GeV^ in the region where data are 
available) in the framework of semiclassical, high-density QCD models and has 
also been found in numerical solutions of the non-linear evolution equations at small x 



13| , p!5| . Let us also point that this model does not correspond to any fixed twist, see 
¥T| for a study of its twist structure. 

The extension of this model to the nuclear case can be made in a straightforward 



manner in the Glauber-Gribov approach [^|: ignoring isospin effects which are negli- 



gible at small x where the model will be applied, we will substitute adp{x,r) by 

(^dAix, r)= d% adA{x, r, 6), (7) 



with h the impact parameter of the center of the dipole relative to the center of the 
nucleus and 



exp ( -]^ATA{h)(ydp{x,r] 



(8) 



the total dipole-nucleus cross section for fixed impact parameter, with adp{x,r) given 
by Eqs. (|) and (^. 7^(6) is the nuclear profile function (longitudinal integral of 
the nuclear density, Ta(6) = (iz p^(z,6)) normalized to unity, / d% TA^b) = 1; 
we employ a nuclear density in the form of a 3-parameter Fermi distribution with 



parameters taken from With this normalization we recover the dipole-nucleon 

cross section making a power expansion for small ATA{b)aQ, keeping the first term and 
putting A = 1. Also the centrality (impact parameter) dependence of the structure 
functions can be computed by direct substitution of crd^(x, r, 6), Eq. (||), into Eq. (|]). 
This model implies a new scaling for nucleus of the type of the r-scaling for proton and 
a new saturation scale, which we will discuss in Section ^. These two Eqs. (|^ and (^ 



constitute the central point in our extension of the model of |35] to the nuclear case. 



^In the 4-flavour version, a flavour dependence is introduced in through variable x and some 
deviation of the scahng may also appear through variable e in the photon wave function due to the 
larger charm mass. 
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Using Eqs. (|T|)-(P) we can compute the nucleon and nuclear structure functions and 
the corresponding ratios, which we will compare with experimental data in the next 
Section. But first let us discuss the region of applicability of the model: This should be 
that of small x (due to the use of the model of for the nucleon, to the neglection of 
isospin effects and to the use of (||) which requires a large coherence length, achieved 
at small x, of the photon fiuctuation), and of small and moderate < 20 GeV^. 
Although the r-scaling may suggest that the region of applicability in could be 
wider, we think that a safe extrapolation to higher would require to implement 



DGLAP evolution Also the extrapolation to very small x could imply effects 

of gluon or pomeron fusion like those included in the non-linear evolution equations 
mentioned previously, so it should be taken with care. Nevertheless, numerical studies 
T3| show that the onset of the non-linear effects is quite smooth, becoming large for 



extremely small x. So we consider^ this model as a reasonable approximation for values 
of a; > 10^^ ^ 10^®, which are those relevant for RHIC and LHC Due to the fact that 
for nuclear structure functions the amount of experimental data is much more limited 
than that for nucleons, we will perform the comparison with experimental data for 
X < 0.02. 



3 Comparison with experimental data 

Here we show the results of the model together with available experimental data. In 
Figs. 1-5, the results of the model in the 3-flavour (solid lines) and 4-flavour (dashed 
lines) versions (see the explanations below Eq. (^) are compared with experimental 
data from H, H, ||, ||. In these Figs., R{A/ B) = [BF2a{x,Q^)]/[AF2b{x,Q% 
Even when joined with a line, the results of the model have been computed at the 
same {{x), (Q^)) as the experimental data; in the latter, the inner error bars are the 
statistical errors, while the outer error bars show statistical and systematic errors 
added in quadrature (with overall normalization uncertainties ignored). Considering 
both the absence of any free parameter to go from the nucleon to the nuclear case 
and the simplicity of the model, we find the agreement quite reasonable. Concretely, 
from Figs. 3 and 4 it can be concluded that the A- and Q^-dependences of the data at 
■^In it is seen that the scaUng induced by the non-linear evolution is fulfilled for values of 
y — b^/iNcCts)]^T^{xo/x) greater than 2.2, which for as — 0.2 and xq — 0.01 means x smaller than 
10-7. 
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fixed X are well reproduced, while in Figs. 1, 2 and 5 the x-dependence is reasonably 
reproduced, taking into account that these latter Figures contain data with going 
from ~ 0.15 GeV^ for the smallest x values to ~ 4.5 GeV^ for the highest x. 

As the version with 3 flavours is simpler, gives an equally reasonable agreement 
with nuclear data than the 4-flavour version, and produces a better description of 



the nucleon data |3§], from now on we will restrict our computations to the 3- flavour 
version. 

Let us turn now to the behaviour of the ratio of longitudinal to transverse cross 
sections. Experimentally |^ large nuclear effects have been observed: the ratio o"l/o"t 



in N (^He) over oi^joq^ in D has been found to reach values as high as ~ 5 (~ 2) for 
0.01 < X < 0.03 and < \ GeV^. Some explanations point to nuclear enhanced 
power corrections, but the experimental data are under reanalysis and the evidence 



of such strong nuclear effects is now dubious. While this important point has to be 
settled, it is clear than in our model such strong effects are not present, as the nuclear 
effects are contained in the Glauber-Gribov cross section in Eqs. (|^) and ® which is 
common to both longitudinal and transverse cross sections, see Eq. (^). In Fig. 6 we 
show the results in our model for the ratio ctl/ctt in nucleus over o"l/o"t in proton, for 
C and Pb. It can be seen that the nuclear effects never go beyond ±12 %, which is a 
clear prediction of our model. 

In Fig. 7 we present predictions of the model for the ratio -F2a/(^-^2p) for C and 
Pb, together with the x-evolution for Be, C, Al, Ca, Fe, Sn and Pb at fixed = 2.25 
GeV^. A clear evolution with can be seen, which in this model is due to the interplay 
between the (transversal and longitudinal) probabilities to get a dipole of size r, Eqs. 
(§) and (H), and the dipole-target cross section, Eqs. (|]), (|^) and (|), and cannot be 



addressed to any concrete twist but to an admixture of all twists, see At large 
enough Q"^ this dependence on will eventually disappear, as this model, as stated 
previously, leads to Bjorken scaling. 

As a last point in this Section, let us comment on other possible options to get 
nuclear structure functions in the framework of the dipole model. A simple form for 
the dipole-nucleus cross section is suggested by high-density QCD |TD|, |I^ |r^: 



o^dA{x,r) = d b 



1 — exp 



(9) 
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We have tried several relations between (the saturation scale in nuclei) and that 
in proton, Q^. On the one hand, we have used a relation coming from the running of 



the coupling, of the type 



On the other hand, we have imposed the first scattering approximation (vahd for r —>■ 0) 
in Eq. (i, 

Ql^ = UTAib)aoQl (11) 

(in this expression the value of the running coupling evaluated at the appropriate scale 
is hidden in ao, see e.g. [|l3l)Q- But so far we have not succeeded in getting a satisfactory 
description of experimental data, not even on a qualitative level: either too strong a 
shadowing is observed or too fast an evolution in x (and too slow in Q"^) is obtained. 
Indeed Eqs. (^ and @ contain different physical assumptions on the nature of the 
scattering centers: while (||) considers multiple scattering on single nucleons (described 
by the saturating form (^), (^ implies scattering on a black area filled with partons 
coming from many nucleons. Our lack of success in reproducing the experimental data 
with (I) suggests that higher order rescatterings are actually needed in the exponent 
of Eq. (^) for the proton, and that the asymptotic region where Eq. @ should be 
valid to describe data on F2A integrated over impact parameter, is not reached yet (i.e. 
the grey region is still dominating the scattering); a very important test for the form 
(^) would be its ability to describe experimental data on diffraction (which is indeed 
fulfilled by (|^) ^5[). In the next section we will address the behaviour of the saturation 
scale in our model. 



4 Unintegrated gluon distribution and saturation 
scale 

As stated in the Introduction, in the fcr-factorization scheme a key ingredient is 
the unintegrated gluon distribution of the hadron, V9^(x, k, h) (sometimes it appears in 
^The fact that both and may have roughly the same x-dependence can be justified by the 
following qualitative argument: is related with p^-broadening in nucleus, — nA{h)Ql, with 
nA{b) the number of scatterings at impact parameter b. As nA{h) oc ATA{h)xGas{Q^gA) I Q^sA (^^^^ ^ 
perturbative QCD cross section evaluated at scale Q^^a ^he gluon distribution in a nucleon) 

and Ql oc xG, both QIa ^^'^ Qs show the same x-behaviour (modulo the logarithm coming from the 
rmining coupling). Special thanks are given to D. E. Kharzeev for discussions and suggestions on all 
these points. 
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the literature as / = k'^({> ||2^, |3T[|), with k the transverse momentum. This ^a{Xi k, b) 
at fixed impact parameter b is related, at lowest order in /cT-factorization to the 
dipole-nucleus cross section by a Bessel- Fourier transform (see []T3|, |TB|): 

V?A(a:, fc, b) = - / fc^ / — exp (ifc ■ r)adA{x, r, b), (12) 

with k"^ = k ■ k, r"^ = f ■ r and vectors defined in the two-dimensional transverse space. 
The unintegrated gluon can be related to the 'ordinary' gluon density (that used in 
coUinear factorization |25||) by 

xG{x,Q'^,b)= dk"^ <^A{x,k,b), (13) 

although this expression must be considered with great it is only true for 

large Q"^ ^ Q'j. (the actual relation is not with the coUinear glue but with the gluon 
distribution in the light-cone wave function of the hadron, see |jlO|, |12|, p!7|). 

For the proton, Eq. ([T2|) leads to the result (pp{x,k) cx cro(A;^/Qs) exp (— fc^/Q^) 
P5[|p|. For the nucleus, using the technique outlined in the Appendix of (or sim- 
ply applying [Nc/ {2Tc'^as)]k'^'Vl to function 0o defined in Eqs. (31) and (34) in that 
Reference), we get 

^A{x,k,b) = — ^^^^^--^exp -— , 14 

with 

B = ^ATAib)ao. (15) 

As in the case of proton, Eq. ( |T4|) shows explicit scaling in k'^/Ql\ besides, in this 
Equation (as in (|^)) the result for proton is recovered making a power expansion for 
small B, keeping the first term and putting there A = 1. For realistic values of 5 < 3, 
Eq. (|l^ turns out to be very suitable for numerical computations, as the convergence 
of the series in n is very fast and only a few terms are needed to get ^^^(a;, k, b) to the 
desired accuracy. 

In Fig. 8 we show the unintegrated gluon distribution for proton, and for Pb in 
three cases: central (6 = 0), peripheral (6 = 7 fm), and integrated over 6, and for two 
values of X = 10~^ and 10~^ (in these computations there is no Q^, so the substitution 
^In some proposals it is considered that the unintegrated gluon distribution should tend 

to a constant as fc — > 0. For discussions on the 'correct' definition and behaviour of this quantity, see 
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in Eq. (|^) is meaningless here and we will make no distinction in this Section between 
X and x\ this also avoids the complication of flavour dependence in case the 4-flavour 
version is to be used, although as stated previously we will use the 3-flavour version). 
A scaling in k'^/Ql^ (with identified with the position of the maximum, see below) 
is perfectly visible, as in each case the curves move to the right with decreasing x while 
keeping their shape and size (this phenomenon has also been found in the framework 
of the non-linear equations for small x |1T0| , |13|, |l5l and the solution called a soliton 
wave [|1^). Besides it can be seen that the shape of the curves for different cases is 
quite close, the only differences being the height, the position of the maximum which 
at fixed x varies from left to right for proton and Pb with increasing centrality, and 
the logarithmic width which slightly increases with increasing centrality (being e.g. at 
1/10 the maximum height, 2.20 for proton and 2.40 for central Pb). 

Keeping in mind the difficulties to identify at small and moderate the integral 
of the unintegrated gluon distribution with the ordinary gluon density, Eq. ([T3|) (see 
the comments below it), it is still tempting to use this Equation and compare with 
other approaches. A comparison at = 5 GeV^ of the results of our model with 
others, for the ratio of gluon densities in Pb over proton, can be found in |2^. There 
it can be seen that our results at x ~ 10^^ roughly coincide with those of [1^, |5^ 



and are higher than those of ||22|, while at x ~ 10^^ they become smaller than 
those of [|19], |5^ , get close to those of and are still larger than those of . Apart 
from the constraints coming from existing DIS experimental data on nuclei which are 



very loose for the glue at small x, in |19] the saturation of gluon shadowing comes 
mainly from the initial condition for DGLAP evolution where this saturation has been 



imposed, while p^ , pO|] are theoretical models and in |51I] the behaviour of the glue has 
been fixed in order to reproduce charged particle multiplicities in AuAu collisions at 
RHIC. Additional caution has to be taken to compare our results with those coming 
from DGLAP analysis ]19[: our ratios for the glue at some moderate, fixed and 
very small x result smaller than the ratios for F2 at the same x, Q^, which leads to 
problems with leading-twist DGLAP evolution as our model mimics the DGLAP 
increase of F2 along a wide Q^-region, this is most probably related with the use of Eq. 
(p^) at too small Q^. In any case, our model leads to Bjorken scaling at some which 
increases with decreasing x: the Q^-independence of the gluon density is achieved in 
our model at ^ 15 GeV^ for x = lO'^ and at ~ 300 GeV^ for x = 10"^. Of 
course this is due to the exponential decrease of the unintegrated gluon at large k"^, 
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see (p!^); other proposals (see e.g. [|T2], |2^) consider a behaviour oc at large fc^ 
which obviously leads to a logarithmic, DGLAP-like increase of the gluon density with 
increasing Q^. 

Now we turn to the saturation scale. While for the case of the proton its definition 
is quite clear in coordinate space, for the nucleus a clean definition is better obtained 
in momentum space, where it corresponds to the maximum of the unintegrated gluon 



distribution |]T2|, [T^, |T5|, The results are shown in Fig. 9 (upper plot) for the same 
cases as in Fig. 8. Contrary to naive expectations, the saturation scale between proton 
and central nucleus is not simply proportional to A^^^ but has a prefactor which makes 
the result smaller (turning the expected factor 5.9 for Pb into a factor 2.0). This can 
be understood using analytical estimates. Taking the exponent in Eq. (||) to be 1/4 
(in analogy to the case of the proton), we get 

-1 



Ql (16) 



2ATA{h)a^-l, 
which gives a 

QIa - 2.4Q2 for Pb at 6 = 0. For AT4(6)ao > 1, we get exactly Eq. (|llj), 
recover the expected proportionality Ql^ oc A^/^Ql and find with this asymptotic form 
QIj^ ~ 2.5(5^. As the form (||) gives a nice agreement with data and (|), (p^T]) do not, we 
conclude that preasymptotic contributions (i.e. not corresponding to ATA{h)aQ oo) 
play a dominant role for the description of present available data. 

To perform some comparison with other approaches, let us see the saturation scale 



coming from the numerical solution of the non-linear evolution equation in [|T5] . There a 
form (for the evolution obtained starting from a phenomenologically motivated initial 
distribution) is obtained that can be approximated, for central (6 = 0) Pb by the 
equation 

g2^^^o.37g^p -6.43 + 0.78 In (^^) ; (17) 

here xq = 0.01, a value of as = 0.2 has been used and the result is in GeV^. The coeffi- 
cients in this Equation coincide with those obtained in other numerical and analytical 



studies [|T2| , [T^, Its results for A = 1 and A = 208 are shown in Fig. 9 (lower plot), 
but one should keep in mind the way in which the coefficients in Eq. (|l^) were obtained 
T5[] : They come from a fit to the position of the maximum of the unintegrated gluon 



distributions for different centralities and values of x but when the scaling induced 
by the non-linear evolution has already set in. So, from the first footnote it can be 
concluded that they correspond to x < 10"^ and thus to a much lower value of x than 
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the region where we have actually plotted them in Fig. 9. From the comparison of the 
two plots in Fig. 9 it is clear that the non-linear evolution produces a much steeper 
behaviour compared with the model in this paper (i.e. an exponent 0.78 compared to 
0.288), while the asymptotic A-dependence is not really different (0.37 ~ 1/3). 



5 Conclusions 

A simple model for nuclear structure functions in the region of small x, and small and 
moderate Q^, has been presented. It is a parameter-free extension of the saturation 
model for the nucleon of p5[ in the Glauber-Gribov approach, although in principle 



any saturating model which correctly describes nucleon data, e.g. could be used. 
This simple extension of the model for the nucleus should be valid until non-linear 
evolution effects |10|, |Tl|, |T^, 0, 0, |T5|, |T6[ become of importance, and for not too 



high where implementation of DGLAP evolution should be required. Conservative 
estimates establish the region of validity of the model to be 10~^ -v- 10~^ < x < 10^^ 
and g2 < 20 GeV^. 

A reasonable agreement with experimental data is obtained for the x-, A- and Q^- 
dependence. The longitudinal-to-transverse cross section ratios show nuclear effects 
smaller than ±12 %. Predictions of the model for values of x smaller than those 
available to present experiments have been given. 

The unintegrated gluon distribution coming from our dipole-nucleus cross section 
has been obtained and studied. Besides the behaviour of the saturation scale in this 
model has been shown and discussed. It turns out to behave oc A^^'^ for large A as 
expected, but large prefactors make the resulting scale for nuclei smaller than naive 
estimates, e.g. for central Pb 2 times bigger, instead of 6 times, than the corresponding 
scale for proton. 

Our model could be used to provide the starting condition for DGLAP evolution, 
as performed in the approaches of |T^, for some initial scale Ql ^ ^qcd- The uninte- 



grated gluon distribution could be employed to compute particle production, using the 



/cy-factorization scheme q^, in high energy collisions involving nuclei p8|, pO 
Work along these directions is in progress. 

As a last comment, our study implies the existence of a scaling for nuclei of the 
same type as that for proton [Q, with a non-trivial relation for the A-dependence 
of the saturation scales between both cases, while the x-dependence turns out to be 
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the same. An experimental, model-independent extraction of the saturation scale in 
nuclei would be very useful to settle the discussions on the relevance of the semiclassical 
approach for existing or future experiments, and the region of validity of perturbative 
QCD and DGLAP evolution. Such an issue would be best explored, and our model 



tested, in high energy lepton-nucleus colliders [p3|. 
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Figure captions: 

Fig. 1. Comparison of the results of the model (see text) in the 3-flavour (solid lines) 
and 4-flavour (dashed line) versions with experimental data at small x for the ratios 
C, Ca and Pb over D, from Q. 



Fig. 2. Id. to Fig. 1 but with data for Be, Al, Ca, Fe, Sn and Pb over C, from |144[ . 
Fig. 3. Comparison of the A-evolution of the results of the model in the 3-flavour 
(open circles) and 4-flavour (open triangles) versions with experimental data at small 



fixed X for Be, Al, Ca, Fe, Sn and Pb over C, from [Q. 

Fig. 4. Comparison of the Q^-evolution of the results of the model in the 3-flavour 
(solid lines) and 4-flavour (dashed line) versions with experimental data for Sn over C 
at small flxed x, from [051. 



Fig. 5. Id. to Fig. 1 but with data for C and Ca over D, from [46 . 
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Fig. 6. Results of the model for the ratio Oi^joT in C and Pb over (Tl/ctt in proton 
versus x. In the plots, lines going from the bottom to the top correspond to = 0.1, 
0.5, 1, 2.25, 5, 10 and 100 GeV^. 

Fig. 7. Results of the model for the x-dependence of F2a/ {AF2p). In the two upper 
plots, results for C (upper plot) and Pb (plot at the middle) versus x are given for 
= 0.1, 0.5, 1, 2.25, 5, 10 and 100 GeV^ (lines going from the bottom to the top). 
In the plot at the bottom, F2a/{AF2p) is drawn versus x at Q"^ = 2.25 GeV^ for Be, C, 
Al, Ca, Fe, Sn and Pb (lines going from the top to the bottom). 

Fig. 8. Upper plot: Results of the model (in GeV"^) for the unintegrated gluon 
distribution for proton (solid line) and for Pb integrated over b (dashed line). Lower 
plot: Results of the model for the unintegrated gluon distribution for peripheral (6 = 7 
fm, solid line) and central (6 = 0, dashed line) Pb. In each case, two curves are provided 
for two values of x = 10^^ (curve to the left) and 10^^ (curve to the right). 
Fig. 9. Upper plot: saturation momentum in the model for proton (solid line), and 
for Pb in three cases: central (6 = 0, dashed-dotted line), peripheral (6 = 7 fm, dashed 
line), and integrated over 6 (dotted line). Lower plot: saturation momentum in the 
numerical solution of the non-linear equation in [^, Eq. (0), for A = 1 (solid line) 
and A = 208 (dashed-dotted line). Notice the difference in vertical scales between the 
plots. 
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